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ABSTMCT 


We have done numerical calculationsof two 
co-located 
nuclear 
bursts 
separatedby 5 s in time. 
The yields were 
both 1.2 Mt. 
At the time of the second burst, the 
central 
density 
in 
the 
fireball has 
dropped 
to 
l/27th of 
atmosphericdensity. 
The 
distance 
to 
which 
energy 
is 
deposited increasesas the cube root of the density so the 
initial radius of the 
hydrodynamic shock 
of 
the 
second 
burst 
is 
about 
three 
times 
that 
of 
the 
first burst. 
Similarly,the reduced density has a profound effect on the 
radiativehistory of the second burst. 
In fact, the second 
burst has a power-timehistory comparable to a single burst 
at 
similar density in the unperturbedatmosphere-- one at 
30-40 km. 
The second burst 
radiates 60% 
of 
Its 
energy 
compared with roughly 28% for the first burst. 


I. 
INTRODUCTION 


We 
have 
done 
a 
numerical simulation of two co-located 1.2-Mt nuclear 


explosionsseparatedby 5 s in time. 
The 
calculations were 
done 
with 
the 


one-dimensional(sphericalsymmetry),coupled radiation-transport/hydrodynamics 


program RADFLO (Zinn 1973). 
We use tables of equation of 
state 
and 
opacity 


data 
for dry air. 
The program has successfullyreproducedatmosphericnuclear 


test data and, with careful attention to numerical resolution,agrees well with 


the 
predictionsof appropriateanalytic results (e.g., jump conditions,blast- 


wave scaling). 


11. 
FIRST BURST 


Single-burstevolutionhas 
been 
described by 
Glasstone 
(1964) , 
Brode 


(1968), Zinn 
(1973) and others. 
For the present set of calculationswe have 


chosen to model a 1000-kg, 1.2-Mt source. 
The 
choice 
is 
rather 
arbitrary. 


Initially, half the energy resides in thermal energy of the device and half as 


kinetic energy of the expanding 
(850 
km/s) 
weapon 
debris. 
Initial energy 


deposition in the ambient air is due to x rays radiated by the hot source. At 


the end of x-ray deposition the 0.3 eV temperaturecontour is at 
100-m 
radius 


although the steepest gradientsof temperatureand pressure are near 30 m. 


‘t’hehot 
core 
grows by radiativeexpanson for the first 100 ps. 
At that 


time the core radius is 56 m and its characteristictemperatureis 60 
eV. 
At 


the 
same 
time the expandingweapon debris has reached 19-m radius. The sound 


speed in the core is 59.4 km/s compared with particle 


debris 
shock; the debris shock is rather weak with a 


sound speed) of 
only 
1.5. 
Profiles of 
velocity, 


temperature at 
100 ps are shown in in Fig. 1. 
Note 


has begun to form at the core edge. 


speeds of 43 km/s in 
the 


Mach number (shock speed/ 


pressure, density, and 


that a hydrodynamicshock 


The hydrodynamicshock forms when 
the 
radiative expansion speed 
drops 


below 
ambient 
sound 
speed. 
Prior 
to 
this 
time 
radiative expansion is 


sufficientlyrapid that no hydrodynamicsignal can keep ahead of the 
radiative 


front. 
At 
100 
ps the core expansion speed is of the order of 0.5 km/s while 


the sound speed just outside the core is 3.5 km/s. 


Formation of the shock marks 
transition to 
the 
hydrodynamic phase 
of 


expansion. 
Characteristics of 
the 
expansionare summarizedin Fig. 2. 
The 


figure includes shock radius (Rs), density (ps), particle velocity 
(Us), 
and 


pressure 
(Ps); 
central 
pressure 
( 
and density (pc); total thermal 
power 


(pTH) 
and the 
integratedradiative 
10SS (-fpT.@t) 
G 


The radiative output remains unaltered by the hydrodynamicshock until the 


shock 
becomes 
optically thick 
at 
about 
0.7 
ms. 
Subsequent brightness 


variationsare determinedby the interplayof shock strength (Zinn and Anderson 


1973), 
fireball area, 
and the optical thicknessof the fireball (Zinn 1973). 


Although the derived power-timecurve during the first maximum is dominated by 


fluctuations of numericalorigin (high-resolutioncalculationsshow very small 


fluctuations),the fireballreaches first maximum near 3.5 ms. 
By 
this 
time 


the fireballhas grown to 133 m and is well into the classicalblast-wavephase 


(Taylor 1950, Brode 1955, Sedov 1959). Examinationof the shock radius history 


in Fig. 2 shows that the fireball entered the blast-wavephase (Rat0“4) at 2 ms 


2 


when 
the 
radius was 110 m, twice the core radius when the hydrodynamicshock 


formed.* 


Profiles of 
the 
state variables (density,pressure, temperature)and of 


velocity are given in Fig. 3 for 40 ms, a time in the midst of 
the 
blast-wave 


phase. 
Although the calculationhas only modest resolutionin the shock front 


(br =60 
cm),jump conditionsare near those expected from 
the 
shock 
velocity 


(vS = 3.34 km/s = 0.4 Rs/c), the ambient sound speed (350 m/s) and the ratio of 


specific heats in the shocked gas (ys = 1.261). For example, the peak pressure 


is low by only 2%. 


Eventually,as the shock wave weakens and the fireball becomes transparent 


(near the time of second maximum = 1 s), fireball evolution is no longer well- 


described by the blast-wavesimilaritysolution. When the shock weakens to the 


point that entrainedair is no longer heated above 0.3 eV, the shock ceases 
to 


be luminous and detaches (“shock breakaway”)from the visible fireball surface. 


This happens at about 40 ms 
in 
the 
present 
calculation. 
Further, as 
the 


central 
pressure drops, 
positive pressure gradients (~ 
> O) create inward 


acceleration,slowing the growth of the visible fireball. The 0.3 
eV 
contour 


begins 
to move inward (albeit at only a few meters per second) at about 1.5 s, 


near the time that the central pressure drops below atmosphericpressure. 
The 


central 
pressure falls to about 75% of ambient at 3 s and then recovers as the 


fireball shrinks. Profiles at 5 
S, 
just before the second burst, are presented 


in Fig. 4. 


111. 
S 
B 


Five 
seconds after the first explosion, an identical source was placed at 


the center of the computingmesh. 
Because expansion of the first fireball had 


*ThiS 
fireball initiallygrew by radiative 
expansion 
and entered 
the blast-wave 
phase when roughly 8 times the mass of the fireball core had been swept 
the 
shock. 
‘5 
by 
Wilke (1982) discusses the evolution of microfireballs(10-1 kt) 
generated by the laser-inducedexplosionsof micron-sizedglass spheres. These 
fireballs, which 
have no radiativeexpansion phase, also begin the blast-wave 
phase when the shock has swept up air about 8 
times 
the 
mass 
of 
the 
glass 
spheres. 


created 
large 
cells 
in 
the 
interior, 
we rezoned 
the 
inner 
mesh to 
give 
cell 


widths 
of 
roughly 
5 m. 


The first burst created an x-ray heated core of roughly 30-m radius, which 


then 
doubled 
In 
size 
by 
radiative expansion. By the time the hydrodynamic 


shock formed at 0.1 ms, the 
core 
radius was 
55 
m. 
We 
would 
expect 
that 


distances 
reached during 
the 
radiativephase would scale inverselywith the 


cube root of the density. With the central density reduced a factor of 27 from 


atmospheric density, we 
would 
expect 
an 
initial core radius of 90 m and a 


radius for hydrodynamicshock formationof 165 m. The temperatureprofiles for 


both 
bursts 
at the end of x-ray depositon and at shock formationare shown in 


Fig. 5. 
The results are close to those expected from simple density scaling. 


Evolution of the second burst begins with the x-ray 
deposition described 


above 
and 
with 
rapid 
expansion of the weapon debris. Figure 6, which gives 


radial contours of density as a function of 
time, 
clearly 
shows 
the 
debris 


shock at roughly 45 m at 0.1 ms after the second burst. 
The hydrodynamicshock 


begins to form at 200 m at about 0.5 ms. 
The debris shock dissipateson a time 


scale 
of 
milliseconds due 
to 
both 
the 
general 
expansion induced by the 


hydrodynamicshock and the infall produced by the rapid depressurizationof the 


center during initial debris expansion. 


The 
density profile inside the pre-existingfireball is fairly flat. 
The 


edge of the density well is several hundred meters from the 
burst 
point. 
In 


such 
a relativelyhomogeneousregion we would expect the second burst to begin 


blast-wave expansion at 
400 m 
-- 
about 
twice 
the 
core 
radius. 
Sedov 


(1959:p.260 ff) 
and 
others have shown that sphericalshockwaves moving into 


regions of variable density (paR‘w) have a shock front radius 
proportional to 
t2.5-w 
In Fig. 6 the shock wave expands as t2/5 
. 
until it encountersthe edge 


of the well, after which it expands more slowly. This is 
entirely consistent 


with the preceding analysis (w < O). 


The 
general 
evolution of 
the 
second 
burst 
is 
summarizedin Fig. 7. 


Perhaps the most outstandingdifferencebetween 
the 
two 
bursts 
is 
the 
far 


greater 
fraction of 
exploslon energy 
radiated away 
by 
the second burst. 


Whereas the first burst radiated 330 
kt 
(27.5%) of 
the 
1200-kt explosion 


energy, 
the 
second 
radiated 720 kt (60%). The two power-timesignaturesare 


quite different. 


4 


The classic power-timecurve of a low-altitudenuclear burst displays two 


peaks 
and 
is 
a 
phenomena of a strong luminous shock in air. 
As the ambient 
. 
density is reduced,a strong shock develops farther and farther away 
from 
the 


burst 
point, 
i.e., 
it 
takes 
longer 
to 
sweep up a few core masses of air. 


Ultimately,for extremely low-densityair, no hydrodynamic shock 
would 
form; 


the 
energy of the burst would go into prompt radiationand thermal energy. 
In 


such conditions the power-timecurve would 
show 
one 
maximum, near 
to, 
and 


decrease after 
that time. 
As the burst height increases then, the power-time 


curve changes from the classic double-peakcharacter to single-peak character. 


This is illustratedin Fig. 8, which shows the calculatedpower-timecurve of a 


medium yield nuclear explosion at 50 km. 
The second burst in this study is 
in 


a 
density 
“ 
env~ronment,which is characteristicof high altitude. The density 


in the well of the first burst at 5 s is similar to atmosphericdensity at 
30 


to 
40 
km. 
The 
power-timecurve for the second burst at 5 s should have the 


character of a burst in this altitude region. Figure 9 is a calculated power- 


time 
curve for a medium yield nuclear burst at 35-km altitude. The similarity 


to the PTH curve of Fig. 7 is striking. Thus, 
with 
less 
energy 
going 
into 


hydrodynamicmotions, more is available to be radiated as thermal power. 
(Note 


the power values in Figs. 8 and 9 


Finally, we address the late state of the double burst. 
How different is 


it 
from 
the 
final 
state 
of 
the 
first 
burst? 
The 
answer 1s -- not very 


different. Figure 10 shows the central 
temperature history 
for 
the 
double 


burst. 
The 
solid 
curve 
rracks 
the 
first 
burst. 
At 5 s the second burst 


produces a sharp spike and subsequentdecay. 
We replot the decay with 
a 
5–s 


time shift as the dashed curve. 


Although 
there 
are 
early 
differences between the curves, these become 


minor after the fireballsbecome transparent (> 1 s). 
The reasons, of 
course, 


are 
that 
both fireballsare rapidly approachingpressure equilibriumwith the 


atmosphere togetherwith the fact that air is a poor 
radiator at 
these 
low, 


late-time temperatures. Similarly the central density returns to a value near 
d 
4 x 10-5 g/cm3 5 
s 
after 
each burst. 
Figure 11 shows profiles of the variables 


5 
s 
after 
the second burst. 
These should be comparedwith Fig. 4. 
Although 


G 


the calculationobviously suffers from poor resolution,particularlyin the two 


shock 
waves, the overall structure (centralvalues, width of the density well, 


location of shock peaks) are expected to be well determined. 


Iv. 
S 


A second, concentricnuclear explosion inside the density well of a 
prior 


explosion behaves normally 
provided 
that 
appropriate density 
scaling 
is 


observed and that the initial expansion occurs far from the walls. 
The 
power- 


time 
history is greatly 


virtually missing, first 


are elevated. The total 
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Fig. 1. 
Flow–field profiles at 0.1 ms for the first 1.2-Mt burst computed 
with RADFLO. The hydrodynamicshockwave has just begun to form 
near 56-m radius. 
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Flow-fieldprofiles of the first burst at 5 S, just before the 
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Pa is ambient pressure. 
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Power-timehistory of a medium yield burst at 50 km. 
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Power-timehistory of the same explosion shown in Fig. 8 but 
calculatedat 35 km. 
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